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Int  roduct ion 

The  solution  of  Neumann's  change-of-phase  problem  (i.e.,  the  solidification 
or  melting  of  a slab,  whose  temperature  is  initially  uniform  and  is  maintained 
constant  at  the  surface  [l])  by  approximate  analytical  means  has  received  some 
attention  in  the  recent  literature  [2j.  Most  of  the  solutions  available,  how- 
ever, present  a single  approximation,  and  are  not  readily  adaptable  to  obtain- 
ing further,  and  hopefully  more  accurate,  approximations.  In  the  present  work 
a method  of  so  doing  is  presented,  in  which  the  multiple  penetration-depth 
technique  of  [3j  is  applied  to  the  formulation  of  Neumann's  problem  in  the 
form  established  in  [4-].  Some  numerical  results  presented  at  the  end  indicate 
that,  at  least  for  the  case  of  constant  properties,  the  proposed  approach  is 
workable  and  satisfactory. 


Analysis 

The  problem  at  hand  refers  to  the  slab>>0,  initially  liquid  at  a uniform 
temperature  (for  the  case  of  solidification),  whose  surface  temperature 

is  maintained  at  the  constant  value  At  any  time  t >o  the 

liquid  occupies  the  region^(t)<>  c «,  , and  the  solid  the  region  0<  * < 

The  temperatures  T and  in  the  liquid  and  solid  respectively  must  satisfy 
the  Fourier  heat-conduction  equation  (with  primes  and  dots  indicating  differen- 
tiation with  respect  to  ? and  t?  respectively): 
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formed  part  of  the  first  author's  M.S.  thesis  in  the  department  of  Civil 
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^0  ABSTRACT  (Continue  on  reeerae  aide  II  neceaamry  end  Identlly  by  block  number) 

An  approximate  method  of  solution  of  Neumann's  change-of-phase  problem 

Is  presented.  The  solution  makes  use  of  a multiple  penetration-deplh 
approach. 
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and  subject  to  the  solid-liquid  interface  conditions 
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with  si(p)=0-  We  note  that  the  solution  is  of  the  form 

A(t)~  1 (5) 

where  "X  is  an  unknown  constant  and  reference  diffusivity. 

Let  now  [4]  two  auxiliary  solut  ionsT  B,  C)  be  defined  as  satisfying 

eq.  (1)  and  the  following  boundary  and  initial  conditions- 

Ts  («.t ; B,T.>  T.  Tt(»,t;T,.cj=C 

T1(y-°;  B-.T«V  T.+  s T[>.o;T,  c).  T7 

L '*•  / ' » 1 (6^ 

B,T«V  Tc+B  Tj-.t-.T.tV  T, 

It  is  then  easily  verified  that  the  solution  of  the  desired  problem  is  obtain* 
by  choosing  the  constants  ^,C  and  Xso  as  to  satisfy  eqs.  (3)  and  (4). 

An  approximate  form  for  the  auxiliary  temperatures  can  he  taken  in  terms  of 
m penetration  depths  ^£t)as  [3]: 

Tfr  fc  8T\  = T+R-  B 7 <L  .-.Yl-  \ 
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and  where  it  is  assumed  that 


(9) 


The  unknowns  l and  are  calculated  by  satisfying  (1)  approximately  by  means 
of  the  "method  of  moments" , i. e . , by  setting 
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For  the  case  of  constant  properties,  it  is  found  [3]  that 

<L . = A.  fy  t l>'° 
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where  the  values  of  A.  and  I)  are  to  be  obtained  numerically  under  the  re- 

V t 

striction  of  eqs.  (9). 

There  now  remain  to  satisfy  the  interface  relations  (3)  and  (4).  The  former 
gives  T _ 


B - 


C = 


and  the  latter  reduces  to  the  form 
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The  solution  of  eqs.  (13)  involves  a trial  and  error  procedure,  in  order  to 
obtain  values  of/i(t)which  lie  (since  for  between  zero  and  the 

largest  of  the^,(t-^.  In  the  cases  of  one  and  two  penetration-depths  , was 
smaller  than  the  smallest  penetration-depth,  so  that  the  step  function  in 
these  equations  could  be  ignored;  but  for  the  case  of  three  penetration-depths 
this  was  not  the  case  and  the  complete  expressions  had  to  be  used.  A numerical 
comparison  of  the  result  obtained  from  one,  two  and  three  penetration  depth 
with  the  exact  values  is  presented  in  the  accompanying  figures,  for  the  case 
in  which  the  solid  and  the  liquid  have  the  same  properties  («<»|),  for  various 
values  of  p . It  is  clear  that  improvements  result  from  the  use  of  several 
penetration  depths,  although  for  this  particular  problem  even  one  penetration 
depth  leads  to  quite  accurate  results. 
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